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Novel vortex structures are found when a thin superconducting film (SC) is covered with a lattice 
of out-of-plane magnetized magnetic dots (MDs). The stray magnetic field of the dots confines 
the vortices to the MD regions, surrounded by antivortices which "crystallize" into regular lattices. 
First and second order transitions are found as magnetic array is made sparser or MD-magnetization 
larger. For sparse MD-arrays fractional vortex-antivortex states are formed, where the crystal- 
symmetry is combined with a non-uniform "charge" distribution. We demonstrate that due to the 
( ant i) vortices and the supercurrents induced by the MDs, the critical current of the sample actually 
increases if exposed to a homogeneous external magnetic field, contrary to conventional SC behavior. 

PACS numbers: 74.78.-w, 74.25.0p, 74.25.Qt, 74.25.Dw. 



The physics of vortex-antivortex pairs in superconduc- 
tors and superfluids has been of general interest for a 
long time. For instance, such pairs are predicted to ex- 
ist in thin superconducting films at finite temperatures 
due to thermal fluctuations Q. Entropy considerations 
show that above the sharply defined Kosterlitz-Thouless 
transition temperature Tkt-, these vortex pairs start to 
unbind, causing the appearance of a finite resistance. Re- 
cently it was found that symmetry-induced antivortices 
can be formed in mesoscopic superconducting polygons 
in a certain parameter-range, such that the vortex- 
antivortex configuration complies with the geometry of 
the polygon. In our recent work, we studied the vortex 
structure of a superconducting film with a single out-of- 
plane magnetized dot on top . The total flux penetrat- 
ing the superconductor equals zero, and vortices cannot 
form in isolation; vortices and antivortices nucleate in 
pairs. A shell vortex structure was observed, with a vor- 
tex nucleus surrounded by an antivortex core (the so- 
called "vortex- molecule" ) with size-magnetization con- 
trolled magic numbers. These vortex configurations re- 
semble the ones of electron dimples on the surface of liq- 
uid helium, electrons in quantum dots, colloidal suspen- 
sions and dusty particles in complex plasmas . 

In the present article, we report further consequences 
of this superconducting Wigner crystallization, in case 
when a regular array of magnetic particles is deposited on 
the superconducting film. Modern advances in microfab- 
rication and characterization techniques 5] have allowed 
an experimental realization of such SC/ferromagnet(FM) 
hybrid systems. Arrays of magnetic particles are po- 
tential devices for applying well-defined local magnetic 
fields, which modulate the order parameter in an underly- 
ing superconductor. Refs. 0, Q (and references therein) 
have explored a plethora of physical effects, including 
matching effects with ordered pinning arrays, where ad- 
ditional pinning contributions arise due to the magnetic 
nature of the pinning centers. 



Here, we investigate the superconducting state of a 
thin SC film with a square array of submicron cubic mag- 
netic dots with perpendicular magnetization (Fig. 1). 
To ensure that MDs and SC are not electronically cou- 
pled, we assume a thin layer of insulating oxide between 
them. We consider cubic MDs, although most of the pre- 
vious experimental work was done on thin FM structures. 
Making the magnetic dots thicker facilitates their mag- 
netizing in the out-of-plane direction and eliminates the 
extreme peak structure in the stray field profile close to 
the dot edge. The general physical behavior of the SC 
drawn out in this article is immune to the MD-thickness. 

The energy difference between the superconducting 
and the normal state, in units of i7^/47r, is 
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where Hq denotes the applied magnetic field (for ex- 
ample, Ho = iimd if magnetic dots are the only field 
source). Eq. ^ is given in dimensionless form, where all 
distances are measured in units of the coherence length 
^, the vector potential A in ch/2e^^ the magnetic field 




FIG. 1: Oblique view of the superconducting film and oxide 
layer (with thicknesses d and /, respectively) underneath a 
regular array (with period W) of cubic magnetic dots. 
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H in Hc2 = cfi/2e^^, and the order parameter ^ in 
y/—a/(3 with a, /3 being the GL coefficients. The min- 
imization of Eq. (QJ leads to weh known GL equa- 
tions which for thin superconductors {d < ^, A) may 
be averaged over the SC thickness. We solve these two 
coupled equations, following a numerical approach pro- 
posed by Schweigert et al. (see Ref. 8j and references 
therein) on a uniform Cartesian grid with typically 10 
points/^ in each direction. In the present case, we took 
for the simulation region a rectangle Wx x W^, where 
Wx = Wy = 16W (i.e. we simulate 16 x 16 supercell, 
see Fig. 1). Periodicity of the SC and the MD-lattice 
is included by applying periodic boundary conditions for 
A and ^ in the form A(r + b^) = A(r) + V7^i(r), and 
^(r + bi) = ^exp(27ri77i(r)/^o) 1^, where hi=x,y are 
the supercell lattice vectors, and rji is the gauge poten- 
tial. These boundary conditions mean that A, ^ are 
invariant under lattice translations combined with spe- 
cific gauge transformations. Since the vector potential 
of a regular array of magnetic dots is periodic by itself, 
we choose r]x = % = 0. If the sample is exposed to 
an additional homogeneous perpendicular magnetic field 
Hext (Ho = Hmd + iiext) wc usc the Landau gauge 
Aext = Hextxey for the external vector potential and 
rjx = HextWxV while r]y = 0. Note that values of Hext 
may not be chosen freely and must fulfill the flux quanti- 
zation per supercell requirement following from the virial 
theorem 

To explore the superconducting state, we start from 
different (randomly generated) initial configurations, in- 
crease/decrease slowly ("sweep up/down") the magneti- 
zation of the MDs M and let the vortex-configuration- 
solution relax to a steady-state one (in principle, it may 
be metastable). For given M, we recalculate the vor- 
tex structure of the film starting from: (i) the previously 
found configuration during the sweep, (ii) Meissner state 
(^ ^ 1) or (iii) the normal state (^ in the whole 
sample) as initial condition. By comparing the energies 
of all found vortex states we determine the ground state 
configuration. The obtained M — W equilibrium vortex 
phase diagram is shown in Fig. 2 for the case of a SC 
film with thickness d = 0.2^ (this is, for example^an ad- 
equate value for a 50nm Pb film at T/Tc = 0.97 0) and 
GL parameter n = 1.2 (approximately corresponding to 
the experimental values found for Pb, Nb, or Al films), 
covered with an oxide layer (thickness I = 0.1^) and an 
array of magnetic cubes with a = D = (see Fig. 1). 

For small distance between MDs, the positive stray 
field under each magnetic dot is compensated by the neg- 
ative fields of the neighboring dots, which decreases the 
amplitude of the magnetic field modulation seen by the 
SC (see Fig. 3). Note that the total flux through the 
SC is always zero. Since the demagnetizing factor of an 
infinite magnetic film is unity, for W = a the magnetic 
field equals zero everywhere and the SC state exists for 
arbitrary value of M. Hence, the magnetization value at 
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FIG. 2: Dependence of the vortex- antivortex configurations 
on the magnetization (M) and period of the magnetic dot 
(MD) lattice (W). Solid lines denote transitions between 
states with different number of vortex-antivortex pairs per 
unit cell (N), while dashed lines indicate second order config- 
urational transformations for fixed N. The vortex structure is 
illustrated by the Cooper-pair density contourplots as insets 
(blue/red - low/high density), where the thin lines outline the 
unit cells of the MD lattice. 



which the S/N transition occurs decreases with increas- 
ing distance between the dots. In this region an unusual 
phenomenon occurs: increased strength (magnetization) 
of the dots drives the SC directly to the normal state, 
although the appearance of vortex-antivortex pairs is ex- 
pected |2|. Namely, even if the critical conditions for 
their nucleation are achieved, there is not enough space 
for stabilizing antivortices in the narrow negative field 
areas. At the same time, the magnetic field under, and 
especially between the MDs is so large (see W = 3^ result 
in Fig. 3) that the SC state is suppressed globally. 

When increasing VF, the compensation effect dimin- 
ishes, the field under MDs increases, suppressing locally 
the SC-state. Between the dots, negative stray field be- 
comes lower (spread over wider area) and superconduc- 
tivity survives. For small magnetization of the dots, the 
magnetic field is only able to suppress the superconduct- 
ing order parameter under the edges of the MDs, where 
the induced currents are maximal Such MDs act 
as pinning centra for external vortices (e.g. resulting 
from the application of an external homogeneous mag- 
netic field). This situation is very similar to the SC 
film perforated with a lattice of antidots. By increas- 
ing the magnetization, the negative fiux between the 
MDs (as well as the positive one under the dots) in- 
creases, and the magnetic field lines can "join" into an 
antivortex (vortex). The critical magnetization and pe- 
riod W for the first nucleated vortex-antivortex pair per 
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FIG. 3: Profile of the MD-magnetic field across the sample 
for different values of the period of the MD-lattice. 

MD (in the ground state) are denoted in the lower part 
of Fig. 2 by a solid line. Using magnetostatic calcu- 
lations, we determined that along this line, the posi- 
tive flux under each MD is approximately constant, and 
equals ^+^^/^o = 1.291 ± 0.012. We found that the 
additional positive magnetic flux necessary for the nu- 
cleation of the next vortex-antivortex pairs is quantized: 
A^+/<^o = 1.07 =b 0.01. These values of threshold fluxes 
weakly depend on the parameters of the SC, but strongly 
on the properties of the stray magnetic field determined 
by the geometrical dimensions of the MDs (i.e. ^q^i 
increases significantly with the size of the dots a). 

For dense lattices, the ant ivort ices are compressed into 
narrow interstitial channels, forced to form regular and 
consequently rigid lattices, where now ant ivort ices are 
"shared" by the neighboring MDs. The rigidness, to- 
gether with the uniform distribution of vortices and an- 
tivortices, makes these ordered vortex structures resem- 
ble ionic crystals. In addition, one finds significant simi- 
larities in the physical mechanisms of crystallization. The 
somewhat simplified theory of cohesion in the ionic (and 
molecular) crystals assumes that the cohesive energy is 
entirely given by the potential energy of classical par- 
ticles localized at equilibrium positions. Because the 
particles in ionic crystals are electrically charged ions, 
the main term in the interaction energy is the interi- 
onic Coulomb interaction. The other contribution comes 
from the strong short-range core-core repulsion due to 
the Pauli principle, without which the crystal would col- 
lapse. Analogy with our system follows from the present 
relations - namely, vortices and antivortices interact anal- 
ogously to ions, except for the absent core-core repul- 
sion, necessary for crystallization. This stabilizing factor 
is brought in our sample by the presence of magnetic 
dots, which effectively keep the vortices and surround- 
ing antivortices apart. Therefore, each antivortex inter- 
acts with a magnetic dot coupled with the vortex un- 
derneath through a Lennard- Jones-like potential, form- 



ing superconducting ion-pairs^ which when brought closer 
together form a two-dimensional ionic crystal. This fasci- 
nating parallel is best illustrated by the vortex-antivortex 

= 1 lattice from Fig. 2, which corresponds to the ion- 
configuration on the surface of a NaCl crystal (and many 
other salts and oxides, e.g. AgBr, PbS, FeO, etc.). 

In A" = 2 crystal (each dot creates a double vortex and 
two antivortices), antivortex dimers are shared between 
the neighboring dots in such a way that each MD is sur- 
rounded by 4 antivortices arranged in a cross. For small 
the adjacent crosses are tilted with respect to each 
other. The tilt angle changes with W (or magnetization 
M, see dashed lines in Fig. 2) and the configuration 
transforms through a second- order phase transition to 
square symmetry (tilt angle zero). This bipartite crystal 
now consists of two sublattices (of vortices and antivor- 
tices), where sites belonging to one lattice are connected 
only to the sites of another (see, e.g., the surface of the 
ReOs crystal). In A^ = 3 lattice, the orientational degree 
of freedom is lost, since antivortices crystallize in a per- 
fect square lattice. With increasing distance between the 
MDs, the crystallization mechanism becomes more influ- 
enced by the inter- ( ant i) vortex interaction than the im- 
posed symmetry, and the square lattice gradually trans- 
forms into a hexagonal one. 

However, if the MDs are set further apart, the an- 
tivortex lattice bonds gradually break, leading to ori- 
ented clusters (rather than a crystal) of vortex-antivortex 
molecules around each MD, like in the case of a sin- 
gle magnetic dot on top of a SC j^. For example, in 
A" = 1 state, the antivortices are no longer in the cen- 
tral interstitial position, but are bound to a particular 
MD. Their relative position is such that it maximizes the 
distance between them. The A^ = 3 molecules have a 
specific orientational order in which molecules under ad- 
jacent columns of MDs are rotated by 60°, due to the 
repulsion of neighboring antivortex-trimers. Note that 
the distribution of antivortices in this case does not obey 
the symmetry of the magnetic potential (cubic MDs). 
This leads to an interesting scenario, where some of the 
vortex-antivortex pairs may annihilate in order to pre- 
serve the square symmetry of the vortex state at the ex- 
pense of the energy, allowing the vortex configuration 
to crystallize again. In such manner, fractional states 
are formed, where some MDs "share" a vortex-antivortex 
pair (the number of pairs per dot A" becomes a rational 
number). In Fig. 4 we show the Cooper-pair density 
plots of two typical fractional states for the square mag- 
netic lattice N = 2-^/2 and A' = 2^/4. Here two species of 
vortices are present (doubly and triply charged), causing 
the adequate rearrangement of singly charged intersti- 
tial antivortices. Although with higher energy than the 
ground-state (for example, AE ^ 278/cTc for 50nm Nb 
or Pb films at temperatures far from Tc), these states are 
metastable and experimentally observable. 

The predicted new vortex configurations can be ob- 



4 



N=27. N=2V, 



V 1 




1 


r ' ' ■ 




- • 


« 


• 






■T : : 

1 






• 


> •—i 

• 


• 


-H 




1^ 




% 
i 


< 

• 


► * 




> 

• 


> H 

• 


• 


< 


► 

• 


— < 




4 


^ 

• 


y- ■*■ i 

:* • * 


f * i 


* * i 

• 


• * 


• 


t 


• 











-20 -15 -10 -5 5 10 15 20-20 -15 -10 -5 5 10 15 20 

x/^ x/^ 

FIG. 4: The |^|^ contourplots of typical fractional vortex- 
ant ivort ex crystals for the 4x4 unit cell simulation region. 

served experimentally by using e.g. scanning probe tech- 
niques like Hall and Magnetic Force Microscopy. These 
vortex-antivortex structures will strongly influence the 
pinning properties and the SC phase diagram. If one 
adds a homogeneous external magnetic field, such that 
the number of additional flux quanta matches the num- 
ber of antivortices at the interstitial sites, one expects 
that annihilation occurs, resulting in a well known match- 
ing vortex configuration with all vortices pinned by the 
MDs. This leads to a peak in the critical current, as a 
function of the applied external field. This mechanism 
explains the recently observed phenomenon of magnetic- 
field-induced superconductivity 10]. In order to ver- 
ify this, we exposed our sample with W/(, = 6.25 to a 
homogeneous magnetic field corresponding to the first 
matching field (one vortex per unit cell) and changed 
gradually the magnetization M, starting each time from 
the normal state. Then we apply current in the x- 
direction as Acx = const, (now Aq = Amd + -^ext + Ac) 
which does not interfere with our boundary conditions. 
When the critical value of Acx is exceeded, the motion of 
(anti)vortices can no longer be prevented and supercon- 
ductivity is destroyed. The results of our calculations for 
the critical current jc as a function of the magnetization 
of the MDs are shown in Fig. 5 for the case with and 
without applied first matching field. If no external field 
is present (dark dots in Fig. 5), higher magnetization M 
induces larger screening currents and jc monotonously 
decreases. The appearance of vortex-antivortex pairs de- 
creases the total current in the sample through the phase- 
change contribution to the current. This leads to some- 
what enhanced critical current which decreases further 
with magnetization and tends to zero. On the contrary, 
in the case of the first matching field (open dots in Fig. 5) 
the critical current equals zero if no pinning is present. 
With increasing magnetization, the antivortex-like cur- 
rents are increased, compensating the current of ex- 
ternal vortices pinned by the dots. For M/Hc2 = 0.418 
maximal compensation is reached, resulting in the max- 
imal critical current. With further increased magneti- 
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FIG. 5: Critical current (in units of jo = cifc2^/47rA^) ver- 
sus the magnetization of the magnetic dot lattice (a/^ = 2, 
W/^ = 6.25) in the case of no applied external magnetic field 
Hext (dark dots) and for the first matching field (open dots). 

zation, the qualitative behavior of jc is similar to the 
Hext = case. Nevertheless, if an external magnetic field 
is present, the critical current of the sample for given 
M is found to be actually higher. This demonstrates 
that, contrary to conventional superconductors, the su- 
perconductivity in SC-FM heterostructures is effectively 
enhanced by an applied magnetic field. 
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Note added in proof: Very recently, Priour and Fertig 
[Phys. Rev. Lett. 93, 057003 (2004)] studied a similar 
system where an extremely thin SC and MDs were in 
immediate vicinity. This leads to much sharper magnetic 
field profiles, enhancing disorder in the mixed state. 
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